The vascular properties of the BOLD signal

The work presented in this thesis is intended to contribute towards the understanding of the cerebral vascular behaviour in response to changes in neuronal activation. The blood oxygenation-level dependent (BOLD) functional magnetic resonance imaging (fMRI) signal provides an indirect measure of neuronal activation, arising from a combination of metabolic and vascular (blood flow and blood volume) changes local to the activation. Therefore the BOLD signal is dependent on local vascular properties as well as on the neuronal activation, leading to a variability of the BOLD signal, based on the underlying vascular structure. It has become an important goal to improve understanding of the mechanisms underlying the BOLD signal in order to separate out this vascular variability from the underlying correspondence with the neuronal activation. The effect of field strength on the temporal characteristics of the BOLD haemodynamic response function is investigated. An earlier BOLD response onset has been measured with increasing static magnetic field strength, consistent with an earlier microvascular (compared with macrovascular) signal response. This result can be used to improve haemodynamic models of the BOLD signal. Hypercapnia, a vasodilator, has been used both to assess the relationship between transverse relaxation and blood oxygenation at 3 and 7 Tesla and to identify vascular heterogeneity between two equivalent brain regions. A tight, linear relationship was found between the level of hypercapnia and transverse relaxation at both 3 and 7 Tesla, whilst the change in transverse relaxation due to hypercapnia increased 2.1 ± 0.5 fold from 3 to 7 Tesla, indicating an approximately linear relationship across field strength. In a separate experiment, a vascular asymmetry was found between the left and right precentral gyri using hypercapnia. This result highlights the need to account for the vascular contribution to the BOLD signal before using this BOLD signal to make comparisons of neuronal activity across brain regions. Finally, an improved model for calibrated BOLD has been proposed and implemented, which requires fewer assumptions than existing approaches. This uses the BOLD response to some task, repeated both at normoxia and hyperoxia. To assess the validity of this model, the effects of paramagnetic oxygen molecules are considered, both dissolved in blood plasma and in airspaces adjacent to the brain. These effects were found to be small, except for in the frontal cortex

Most small cerebral cortical veins demonstrate significant flow pulsatility: a human phase contrast MRI study at 7T

Phase contrast MRI has been used to investigate flow pulsatility in cerebral arteries, larger cerebral veins and the cerebrospinal fluid. Such measurements of intracranial pulsatility and compliance are beginning to inform understanding of the pathophysiology of conditions including normal pressure hydrocephalus, multiple sclerosis and dementias. We demonstrate the presence of flow pulsatility in small cerebral cortical veins, for the first time using phase contrast MRI at 7 Tesla, with the aim of improving our understanding of the haemodynamics of this little-studied vascular compartment. A method for establishing where venous flow is pulsatile is introduced, revealing significant pulsatility in 116 out of 146 veins, across 8 healthy participants, assessed in parietal and frontal regions. Distributions of pulsatility index and pulse waveform delay were characterized, indicating a small, but statistically significant (p<0.05), delay of 59±41 ms in cortical veins with respect to the superior sagittal sinus, but no differences between veins draining different arterial supply territories. Measurements of pulsatility in smaller cortical veins, a hitherto unstudied compartment closer to the capillary bed, could lead to a better understanding of intracranial compliance and cerebrovascular (patho)physiology

Vascular physiology drives functional brain networks

We present the first evidence for vascular regulation driving fMRI signals in specific functional brain networks. Using concurrent neuronal and vascular stimuli, we collected 30 BOLD fMRI datasets in 10 healthy individuals: a working memory task, flashing checkerboard stimulus, and CO2 inhalation challenge were delivered in concurrent but orthogonal paradigms. The resulting imaging data were averaged together and decomposed using independent component analysis, and three “neuronal networks” were identified as demonstrating maximum temporal correlation with the neuronal stimulus paradigms: Default Mode Network, Task Positive Network, and Visual Network. For each of these, we observed a second network component with high spatial overlap. Using dual regression in the original 30 datasets, we extracted the time-series associated with these network pairs and calculated the percent of variance explained by the neuronal or vascular stimuli using a normalized R2 parameter. In each pairing, one network was dominated by the appropriate neuronal stimulus, and the other was dominated by the vascular stimulus as represented by the end-tidal CO2 time-series recorded in each scan. We acquired a second dataset in 8 of the original participants, where no CO2 challenge was delivered and CO2 levels fluctuated naturally with breathing variations. Although splitting of functional networks was not robust in these data, performing dual regression with the network maps from the original analysis in this new dataset successfully replicated our observations. Thus, in addition to responding to localized metabolic changes, the brain’s vasculature may be regulated in a coordinated manner that mimics (and potentially supports) specific functional brain networks. Multi-modal imaging and advances in fMRI acquisition and analysis could facilitate further study of the dual nature of functional brain networks. It will be critical to understand network-specific vascular function, and the behavior of a coupled vascular-neural network, in future studies of brain pathology

Arterial CO2 fluctuations modulate neuronal rhythmicity: Implications for MEG and fMRI studies of resting-state networks

A fast emerging technique for studying human resting state networks (RSNs) is based on spontaneous temporal fluctuations in neuronal oscillatory power, as measured by magnetoencephalography. However, it has been demonstrated recently that this power is sensitive to modulations in arterial CO2 concentration. Arterial CO2 can be modulated by natural fluctuations in breathing pattern, as might typically occur during the acquisition of an RSN experiment. Here, we demonstrate for the first time the fine-scale dependence of neuronal oscillatory power on arterial CO2 concentration, showing that reductions in alpha, beta, and gamma power are observed with even very mild levels of hypercapnia (increased arterial CO2). We use a graded hypercapnia paradigm and participant feedback to rule out a sensory cause, suggesting a predominantly physiological origin. Furthermore, we demonstrate that natural fluctuations in arterial CO2, without administration of inspired CO2, are of a sufficient level to influence neuronal oscillatory power significantly in the delta-, alpha-, beta-, and gamma-frequency bands. A more thorough understanding of the relationship between physiological factors and cortical rhythmicity is required. In light of these findings, existing results, paradigms, and analysis techniques for the study of resting-state brain data should be revisited

Changes in arterial cerebral blood volume during lower body negative pressure measured with MRI

Cerebral Autoregulation (CA), defined as the ability of the cerebral vasculature to maintain stable levels of blood flow despite changes in systemic blood pressure, is a critical factor in neurophysiological health. Magnetic resonance imaging (MRI) is a powerful technique for investigating cerebrovascular function, offering high spatial resolution and wide fields of view (FOV), yet it is relatively underutilized as a tool for assessment of CA. The aim of this study was to demonstrate the potential of using MRI to measure changes in cerebrovascular resistance in response to lower body negative pressure (LBNP). A Pulsed Arterial Spin Labeling (PASL) approach with short inversion times (TI) was used to estimate cerebral arterial blood volume (CBVa) in eight healthy subjects at baseline and -40mmHg LBNP. We estimated group mean CBVa values of 3.13±1.00 and 2.70±0.38 for baseline and lbnp respectively, which were the result of a differential change in CBVa during -40 mmHg LBNP that was dependent on baseline CBVa. These data suggest that the PASL CBVa estimates are sensitive to the complex cerebrovascular response that occurs during the moderate orthostatic challenge delivered by LBNP, which we speculatively propose may involve differential changes in vascular tone within different segments of the arterial vasculature. These novel data provide invaluable insight into the mechanisms that regulate perfusion of the brain, and establishes the use of MRI as a tool for studying CA in more detail

Global intravascular and local hyperoxia contrast phase-based blood oxygenation measurements

AbstractThe measurement of venous cerebral blood oxygenation (Yv) has potential applications in the study of patient groups where oxygen extraction and/or metabolism are compromised. It is also useful for fMRI studies to assess the stimulus-induced changes in Yv, particularly since basal Yv partially accounts for inter-subject variation in the haemodynamic response to a stimulus. A range of MRI-based methods of measuring Yv have been developed recently. Here, we use a method based on the change in phase in the MR image arising from the field perturbation caused by deoxygenated haemoglobin in veins. We build on the existing phase based approach (Method I), where Yv is measured in a large vein (such as the superior sagittal sinus) based on the field shift inside the vein with assumptions as to the vein's shape and orientation. We demonstrate two novel modifications which address limitations of this method. The first modification (Method II), maps the actual form of the vein, rather than assume a given shape and orientation. The second modification (Method III) uses the intra and perivascular phase change in response to a known change in Yv on hyperoxia to measure normoxic Yv in smaller veins. Method III can be applied to veins whose shape, size and orientation are not accurately known, thus allowing more localised measures of venous oxygenation. Results demonstrate that the use of an overly fine spatial filter caused an overestimation in Yv for Method I, whilst the measurement of Yv using Method II was less sensitive to this bias, giving Yv=0.62±0.03. Method III was applied to mapping of Yv in local veins across the brain, yielding a distribution of values with a mode of Yv=0.661±0.008

The effect of isocapnic hyperoxia on neurophysiology as measured with MRI and MEG

The physiological effect of hyperoxia has been poorly characterised, with studies reporting conflicting results on the role of hyperoxia as a vasoconstrictor. It is not clear whether hyperoxia is the primary contributor to vasoconstriction or whether induced changes in CO2 that commonly accompany hyperoxia are a factor. As calibrated BOLD fMRI based on hyperoxia becomes more widely used, it is essential to understand the effects of oxygen on resting cerebral physiology. This study used a RespirActTM system to deliver a repeatable isocapnic hyperoxia stimulus to investigate the independent effect of O2 on cerebral physiology, removing any potential confounds related to altered CO2. T1-independent Phase Contrast MRI was used to demonstrate that isocapnic hyperoxia has no significant effect on carotid blood flow (normoxia 201 ± 11 ml/min, -0.3 ± 0.8 % change during hyperoxia, p = 0.8), whilst Look Locker ASL was used to demonstrate that there is no significant change in arterial cerebral blood volume (normoxia 1.3 ± 0.4 %, -0.5 ± 5 % change during hyperoxia). These are in contrast to significant changes in blood flow observed for hypercapnia (6.8 ± 1.5 %/mmHg CO2). In addition, magnetoencephalography provided a method to monitor the effect of isocapnic hyperoxia on neuronal oscillatory power. In response to hyperoxia, a significant focal decrease in oscillatory power was observed across the alpha, beta and low gamma bands in the occipital lobe, compared to a more global significant decrease on hypercapnia. This work suggests that isocapnic hyperoxia provides a more reliable stimulus than hypercapnia for calibrated BOLD, and that previous reports of vasoconstriction during hyperoxia probably reflect the effects of hyperoxia-induced changes in CO2. However, hyperoxia does induce changes in oscillatory power consistent with an increase in vigilance, but these changes are smaller than those observed under hypercapnia. The effect of this change in neural activity on calibrated BOLD using hyperoxia or combined hyperoxia and hypercapnia needs further investigation

Calibrated BOLD using direct measurement of changes in venous oxygenation

Calibration of the BOLD signal is potentially of great value in providing a closer measure of the underlying changes in brain function related to neuronal activity than the BOLD signal alone, but current approaches rely on an assumed relationship between cerebral blood volume (CBV) and cerebral blood flow (CBF). This is poorly characterised in humans and does not reflect the predominantly venous nature of BOLD contrast, whilst this relationship may vary across brain regions and depend on the structure of the local vascular bed. This work demonstrates a new approach to BOLD calibration which does not require an assumption about the relationship between cerebral blood volume and cerebral blood flow. This method involves repeating the same stimulus both at normoxia and hyperoxia, using hyperoxic BOLD contrast to estimate the relative changes in venous blood oxygenation and venous CBV. To do this the effect of hyperoxia on venous blood oxygenation has to be calculated, which requires an estimate of basal oxygen extraction fraction, and this can be estimated from the phase as an alternative to using a literature estimate. Additional measurement of the relative change in CBF, combined with the blood oxygenation change can be used to calculate the relative change in CMRO2 due to the stimulus. CMRO2 changes of 18 ± 8% in response to a motor task were measured without requiring the assumption of a CBV/CBF coupling relationship, and are in agreement with previous approaches

Field strength dependence of grey matter R2* on venous oxygenation

The relationship between venous blood oxygenation and change in transverse relaxation rate (ΔR2 *) plays a key role in calibrated BOLD fMRI. This relationship, defined by the parameter β, has previously been determined using theoretical simulations and experimental measures. However, these earlier studies have been confounded by the change in venous cerebral blood volume (CBV) in response to functional tasks. This study used a double-echo gradient echo EPI scheme in conjunction with a graded isocapnic hyperoxic sequence to assess quantitatively the relationship between the fractional venous blood oxygenation (1-Yv) and transverse relaxation rate of grey matter (ΔR2 * GM), without inducing a change in vCBV. The results demonstrate that the relationship between ΔR2 * and fractional venous oxygenation at all magnet field strengths studied was adequately described by a linear relationship. The gradient of this relationship did not increase monotonically with field strength, which may be attributed to the relative contributions of intravascular and extravascular signals which will vary with both field strength and blood oxygenation

The venous contribution to sodium MRI in the human brain

Purpose: Sodium MRI shows great promise as a marker for cerebral metabolic dysfunction in stroke, brain tumor, and neurodegenerative pathologies. However, cerebral blood vessels, whose volume and function are perturbed in these pathologies, have elevated sodium concentrations relative to surrounding tissue. This study aims to assess whether this fluid compartment could bias measurements of tissue sodium using MRI. Methods: Density‐weighted and B1 corrected sodium MRI of the brain was acquired in 9 healthy participants at 4.7T. Veins were identified using co‐registered 1H T*2‐weighted images and venous partial volume estimates were calculated by down‐sampling the finer spatial resolution venous maps from the T*2‐weighted images to the coarser spatial resolution of the sodium data. Linear regressions of venous partial volume estimates and sodium signal were performed for regions of interest including just gray matter, just white matter, and all brain tissue. Results: Linear regression demonstrated a significant venous sodium contribution above the underlying tissue signal. The apparent venous sodium concentrations derived from regression were 65.8 ± 4.5 mM (all brain tissue), 71.0 ± 7.4 mM (gray matter), and 55.0 ± 4.7 mM (white matter). Conclusion: Although the partial vein linear regression did not yield the expected sodium concentration in blood (~87 mM), likely the result of point spread function smearing, this regression highlights that blood compartments may bias brain tissue sodium signals across neurological conditions where blood volumes may differ